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SUMMARY 

A method for the theoretical determination of the lifetime of com- 
pressed plates at elevated temperatures is presented* In this approach; 
linearized equations are used throughout with the assumption that the 
plate material is a standard linear solid. The critical time (lifetime) 
Is determined by reducing the time-dependent behavior to the time- 
independent response of purely elastic buckling. 

Theoretically predicted lifetimes of 2024-T3 (formerly 24S-T3) 
aluminum- alloy plates at 450° F are compared with experimental values 
obtained in previous work. 


INTRODUCTION 

With the increase of flight speeds; it becomes more and more impor- 
tant to investigate extensively the behavior of compressed load- carrying 
structural components such as columns and plates. At high speeds one 
of the effects of aerodynamic heating consists in changing the properties 
of many structural materials from elastic to viscoelastic. Accordingly; 
numerous papers dealing with the problem of such time-dependent be- 
havior; particularly creep bending and buckling; appeared during the 
last decade. A review of these and related questions was presented re- 
cently by Hoff (ref. l) and Shanley (yef. 2). 

As has been repeatedly observed (refs. 1 and 2 ), one of the essential 
features of such time- dependent response is the existence of a critical 
time at which the structure fails under constant compressive loads. A 
variety of viewpoints may be adopted to determine analytically this criti- 
cal time. For example; expressions for the deflection of a compressed 
element can be found; and the critical time can be defined as the time 
at which this deflection becomes infinite. Since in this case the use 
of a nonlinear stress-strain relation (or creep law) is required and since 
the experimental creep curve for metals is nonlinear; the use of nonlinear 
stress- strain relations possibly may be required in order to determine 
analytically any critical time. 
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Another approach could consist in assuming an initial imperfection 
that would increase in time* The critical time (or lifetime) could then 
be defined; still with a linear stress-strain relation; as the time at 
which a certain allowable quantity (such as stress; strain, or deflection) 
is reached. 

Still another point of view is adopted in the present study. It is 
assumed that the plate is perfectly flat and that the constant compres- 
sive stresses will cause time-dependent decrease of Young 1 s (or shear) 
modulus, which in turn is proportional to the flexural rigidity. This 
decrease is determined by developing the plate equations on the basis of 
stress-strain relations derived from spring- dashpot models. The critical 
time is defined as the instant at which the flexural rigidity reaches a 
value that would make a corresponding elastic plate unstable under the 
given loading. 

This approach resembles somewhat the use of the secant modulus in 
the analysis of plastic buckling and is difficult to justify on purely 
theoretical grounds. The main advantage is that it conceptually reduces 
the time-dependent behavior to a purely elastic, time- independent response 
without introducing any new buckling or failure criteria. Therefore, a 
high degree of idealization is embodied, and the physical reasonableness 
can be tested only by comparison with experimental results; but its merits 
for the aircraft structural engineer are obvious. 

Another feature of the present approach, which contrasts it to the 
other methods, consists in the fact that no initial imperfections have to 
be assumed. This is of some significance, since experimental results 
(ref. 3) show that initial deviations from flatness do not influence 
appreciably the lifetime of compressed plates. 

The basic equations governing the behavior of viscoelastic plates 
in compression are derived and discussed in the first section of this 
report. The second section establishes the expressions relating the 
applied compressive force to the lifetime of the plate. 

In the third and last section a quantitative comparison of the theo- 
retical predictions and experimental results is carried out. The material 
constants necessary for such a comparison were extracted from a completely 
different source (ref. 4) than that giving the experimental results of 
the lifetime of plates (ref. 3). 

Because of the scarcity of information on material properties and 
lifetimes and because of the large scatter to be expected in such experi- 
ments, this first quantitative comparison, which gives the correct order 
of magnitude and the correct dependence on the parameters, may be con- 
sidered as promising. 
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SYMBOLS 


a,b 

De 

E 

F 

*h > F l 
G 


h 


K 


k 


MxjMyj^xy 




t 

^cr 

w 


x,y,z 


V 

^d 

e x> 6 yj e z>Vxy 

T 


edge lengths of plate 

elastic flexural rigidity of plate 

upper and lower flexural rigidity of plate, respectively 
Young’s modulus 

uniform edge thrust per unit length of plate 
elastic critical load 

upper and lower critical load, respectively 
shear modulus 
plate thickness 
bulk modulus 

numerical factor, buckling value 
plate moments defined by equation (2) 
linear differential operators 
time 

lifetime (critical time) 
deflection of plate 
rectangular Cartesian coordinates 
Poisson’s ratio 
strain deviator 
strain components 

coefficient of viscosity introduced in figure 2(b) 
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moduli of rigidity introduced in figure 2(b) 
stress deviator 


M-IjM-2 
a d 

cx x , a yj a z> T xy stress components 

time of relaxation .of stress under constant strain 

O 

^ a time of relaxation of strain under constant stress 

Subscripts : 

cr critical 

e elastic plate 

h upper 

l lower 

r relaxed (modulus) 

u unrelaxed (modulus) 


BASIC EQUATIONS 

A flat y rectangular plate with edge lengths a and b and thickness 
h is referred to the system of coordinates shown In figure 1 and is 
loaded along the edges x = 0 ; a by a uniformly distributed (compressive) 
force per unit length F. All four edges are simply supported. 

The stress equation of equilibrium In the transverse direction z 
is (e.g.j see ref. 5) 


Ac + 2 ^ 


F 2?v 


c)x 2 ’ “ ' Sy 2 ‘ £>x 2 

where the moments are defined in the usual way: 

Mx = 


(1) 


rh/2 

rh /2 

rh/2 

I a x z dz >* ^7 = 

V dz; “xy = 


;-h /2 

l-h/2 

1- h /2 


xy 


( 2 ) 



5 


and w is the deflection. In the present linearized analysis, the 
force F is not related to w. To express equilibrium (eq. (l)) in 
terms of one single unknown (deflection w) , the stress-strain relations 
of the material and the geometrical strain-displacement relations have 
to be specified. 


In establishing the relations between the moments M x , My, M^y 

and the deflection w, assume that the plate material is linearly visco- 
elastic. In the absence of volume viscosity, the stress- strain relations 
for such a material may be written down in the form (see ref. 6): 


Pa, = 2Qe, 
d d 


+ o, r + a„ = 


E 


x w y z l - 2v v x y z 


(e + e + e ) 


(3) 


where a^ is the stress deviator, Is the strain deviator, a x , a , 

and a z are the normal components of stress, and e x , Sy, and e z are 

the elongations. The quantity E/3(l - 2v) is the bulk modulus K, and 
P and Q are the linear differential operators in time: 


P = a 0 + a l Jb + 

Q = b Q + b x ^ + 


+ a 


m 


dt m 


+ b 


n 



(4) 


The stress-strain relations (eq. (3)) involving normal stresses and 
strains may be written out in scalar form as 


Pa 


Pa 


Pa 


x 


y ~ 


z 


E P 2 ^ 

•N 


e + 2Q& V 

1 - 2v 3 3 * 

^ x 

E 1-1*1 

e + 2Qe xr 


y 

E 1 - 4 <»' 

6 + 2Q 6^; 

1 - 2v 3 3 

J 


( 5 ) 
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By eliminating e z in the preceding set and neglecting a z , 



( 6 ) 


The procedure is analogous to the treatment of the elongation e z and 

the stress component a in case of elastic plate materials (e.g., see 
ref. 7). 

By assuming, further, the stress-displacement relations to be those 
of the classical linear plate theory, namely 


£>x 2 ' 


e x ■ - z —5 £ V = 


b^y 

by 


b^ w 


2’ r xy - - z ax dy 


(7) 


and by noting that in shear 


^ T xy ~ ^xy 


the moment-deflection relations can be found by combining equations (2), 
(3) , and ( 7) : 


( E P 4 Q \ M _ j£ l 2E P 2 A S 2 w h 3 / E P 2 A OA d 2 v 
\1 - 2v 3 + 3 3 12 \1 - 2v 3 + 3 y 2Q ^x 2 12 \1 - 2v 3 " 3 y 2Q ^ y 2 

/ E P . 4 A __ h 3 ( 2E 2 A d 2 w h 3 / E P 2 A d 2 v 

\1 - 2v 3 3 y M y " 12 [l - 2v + 3 2Q ^ y 2 " 12 \1 - 2v 3 3 2Q ^ x 2 

r, ( E P .4 A „ 2h 3 / E P 4 A A S 2 w 

P (l - 2v 3 3 ^y = 12 ^1 - 2v 3 + 3 ^ 5F5y 


( 8 ) 


then the equilibrium condition, equation (i) in terms of w only, is 
obtained by substituting relations (8) into equation (l): 


(. — i — 1 + §) 

\1 - 2v 3 3/ 


Qh 3 V 4 w + 3P 


( e p 4 A _ 

V - 2v 3 3 7 Sx 2 " 


0 


( 9 ) 
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To treat any particular initial- value problem, the specific form of 
the operators P and Q has to be prescribed. In case of the so-called 
four-parameter model (fig. 2(a)), which was discussed and applied, for 
example, by Glauz and Lee (ref. 8), these operators are of the form 



where pq, T]q, p 2 j an ^ ^2 are ^e four constants of the model. 


( 10 ) 


Since in the case of a four- parameter model the operators P and 
Q are of second order, considerable complications are introduced in 
solving the equilibrium equation (9) together with the appropriate initial 
conditions. Since, in the process considered presently, no unloading 
takes place, the omission of the element representing unrecoverable flow 
(t]^ in fig. 2(a)) appeared to be justified. The model was thus reduced 

to a three -parameter one, as shown in figure 2(b) . This model was dis- 
cussed extensively by Zener (ref. 9) and is referred to as the one 
characterizing a standard linear solid. 

The operators P and Q are now of first order: 



Several other quantities, defined in terms of the model constants 
p-p P 2 > and T) may be introduced for convenience and better physical 

insight. In particular, it is useful to know the time-dependence of E 
and v in terms of the shear modulus G. The quantity T e = q/pq is 

the time of relaxation of stress under constant strain, T a = T|/pq + r)/p 2 
is the time of relaxation under constant stress, G^ = p 2 is the relaxed 
shear modulus, and G u = Pq + p 2 is the unrelaxed shear modulus. It 
follows that 


G 

G 


u 

r 


( 12 ) 
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By introducing as the unrelaxed Young* s modulus; E r as the relaxed 

Young T s modulus; as the unrelaxed bulk modulus; and K^. as the re- 

laxed bulk modulus; the absence of volume viscosity requires that 

Ku = Kr = K 


But 


_ 1 E(t) G(t) __ 1 ^r ^r _ 1 ^u ^u 

3 3G(t) - E(t) “ 3 3G r - E r " 3 3G U - E u 


Since (ref. 9) 


then 


Evidently 


E u = 2(1 + v) G u 

_3_ _1_ 1 - 2v 

E r " G r " ^ 


E r / 2(1 + v ) G r 


since, from equation (13), 

Eft) ^ G(t) 

E(t) - (1 - 2v ) G(t) + E u 


(14) 


(15) 


( 16 ) 


Some authors (refs. 10 and ll) prefer to use the form 

E(t) - 2 [1 + v(t)] G u (17) 

By comparison with equation (16), v (t) must be 

/iX (2 + 5v) G(t) - E u 
(t) “ (1 - 2v) G(t) + E u 

It may be noted that; while v(0) = v, v(.) / 1/2 except for a Maxwell 
model (pg - 0) . Further; in the expression 

a x + °y + a z " 3K ( 6 x + 6 y + e z) 



(19) 
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2(1 + v)G u 

the value 3K may be conveniently replaced by — > which is a 
constant. The quantity v has here its initial value. 

Even though for a three-parameter model the operators P and Q 
are of first order, equation (9) on the deflection w will be of second 
order. 

An attempt could be made to lower this order while the same basic 
model is still being considered for the stress- strain relations. This, 
in fact, can be done by extending the elastic stress-strain relations for 
a plate directly to the viscoelastic ones without ever considering, as 
was done before, the three-dimensional stress-strain relations. The 
physical significance of this latter procedure, as contrasted to the 
previous one, will be discussed at the end of the next section. 


In the present case of a three-parameter model, therefore, 


( X + ~ ( G + T a + ve y) 

( x + T e It) 0 y = ( 2 + T a |^)( e y + v£ x) 

( X + x £ Jb) T xy = G r ^ + T cr H) 


This leads to the equation 


h 3 2G r 


12 1 - v 


xy 


( X + T a Je) v4w + ( X + T e F = 0 


or 


— ^ — V 4 v + PF 1-^ =* 0 

12 1 - v v w + n Sx 2 


( 20 ) 


( 21 ) 


( 22 ) 


which replaces equation (9) and is of the first order in time. 


CONCEPT OF ELASTIC BUCKLING APPLIED TO VISCOELASTIC PLATES 

The external force F is assumed to act in the form of a step func- 
tion and keeps its constant value after the beginning of loading at a 
time t = 0. 
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One of the features of the linearized elastic stability analysis is 
that no relation is sought between the force F and the deflection w 
(ref. 12 contains a careful and detailed analysis of the stability 
phenomenon) . 

In the case of an elastic plate, the critical force F g is found 
to be (ref. 12) 

D it 2 

F e = -V k (23) 

b^ 


where D g is the flexural rigidity of the elastic plate 

G e h 3 

D e * 6(1 '-TF 


(24) 


and G e is the shear modulus of the elastic plate* k is a numerical 

factor, the buckling value, the magnitude of which depends on the ratio 

a/b . 


In the case of a viscoelastic plate characterized by the three- 
parameter model (fig. 3), the following statement can be made: Since the 

plate is assumed to be perfectly flat, the constant force F produces 
only constant compressive stresses. Further, at the instant of loading 
the viscosity element is as yet ineffective, and the rigidity is governed 
by the unrelaxed modulus G u . On the other hand, if the force F is 

permitted to act for an infinitely long time, the spring p 2 becomes in- 
effective and the rigidity is governed by the relaxed modulus G r . 

Thus an upper critical load F^, under which the plate buckles in- 
stantaneously, and a lower critical load Fj , under which the plate buck- 
les after an infinite time, may be defined. Hence, for any load F > F^ 
the plate buckles instantaneously* for F< Fj the plate never buckles* 
and for F j < F < F^ the plate buckles after a finite period of time, 
which may be referred to as the lifetime or critical time t cr . 

The expressions for F h , Fp and t cr may be determined readily. 

To calculate F h , let P = 1 and Q = G u , and obtain from equation (9) 



( 25 ) 
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where 


D, 


G u h^ 


h 6(l - v ) 

To find F^, let P = 1 and Q = G r , and obtain similarly 

ty 


D nit 2 


(26) 


(27) 


where 


2(1 + v)G u 

G r h 3 

_(1 - 2v)G r X _ 

3 

2(1 + v)G u 

- 4 

- 


(1 - 2v)G r ■* 


(28) 


To calculate the critical time t cr associated with a load F Z < F < F h , 

proceed as follows: Since the (constant) force F produces constant 
stresses in the perfectly flat plate, the stress-strain relation (eqs. 
(20)) takes the form 


x 


(L 



(29) 


During the interval of time 0 < t < the shear modulus G(t) will 

decrease from G u to G r . Since, at t « 0, t = G u y, then from equation 

(29) 


G(t) _ 

G u G u t/x a Gr 
G 


(30) 


The load F as a function of critical time is therefore (again with the 
use of eq. (9) ) : 


F ^ t cr) _ Dft) 
F h " Dh 


2(1 + v) + (1 - 2 v) 
~u 

2(1 + v) + 4(1- 2v) —P 


2(1 - 


v) 



( 31 ) 
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A plot of G(t)/G u against a dimensionless time t/x a is presented 

in figure 3. Poisson's ratio was taken as v = 0.3. One curve is for 
G r /G u = 0.9 (low viscosity) and the other for G r /G u = 0.1 (high 

viscosity). Figure 4 shows a plot of F(t cr )/F h against t/x a , again 
for G r /G u — 0.9 and G r /G u * 0.1. 


With the simpler equation (22), the same relation (30) is obtained 
for G(t) /G u ; but a different one is obtained for F(t cr )/F h , namely 


F(t) _ D(t)_ _ G( t) 


r h 




(32) 


which now replaces equation (3l). Thus, with the simpler equation (30), 
the curve F(t)/F tl is identical to G(t)/G u , whatever the value of 

a rK- 


It may be of interest to establish the physical meaning of stress- 
strain relations (20) as contrasted to stress-strain relations (3) . A 
consequence of the former is equation (3l), while a consequence of the 
latter is equation (32). Thus, in order to pass from equation (3l) to 
equation (32) , 


2(l + v) + (l - 2 v) 




G,. 


2(1 - v) 


2(1 + v) + 4(1 - 2v) 


= 1 


(33) 


Introducing into this expression the bulk modulus K gives 


K = 


E 


3(1 - 2v ) 


__ 2(l + 

= 3(1 - 


2vJ 


and solving for K gives 


K = 


(1 + V 1 

(1 - 2 


2^7 


G(t) 


(34) 


(35) 


or 


K(il = G(t} 

Ku Gu 


( 36 ) 
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that is, the bulk modulus changes in proportion to the shear modulus. 
Thus, while equations (3) express absence of volume viscosity 
(K = constant), equations (20) describe a material with both volume and 
shear viscosity. 


COMPARISON OF THEORETICAL PREDICTIONS 
AND EXPERIMENTAL RESULTS 

In order to predict a numerical value for the critical time, three 
material constants corresponding to the three elements of the model must 
be known. These constants may be obtained from a creep test. 

Figure 5 shows an approximate reproduction of a curve from reference 
3, which gives the total deformation in percent against time in hours 
for a 2024-T3 aluminum sheet at 450° F subjected to a compressive constant 
stress of 26,000 pounds per square inch. From this plot the unrelaxed 
Young* s modulus E u is readily obtained by dividing the stress by the 
instantaneous strain: 


_ 26,000 « * / 7 7 \ 
E = — - = 6.50x10° psi (37) 

4x10"° 

If the plot is extrapolated by assuming that the total deformation will 
reach 1.6 percent for very large times, the relaxed Young's modulus is 
obtained: 


E = 


26,000 

1.6x10“ 


1.60x10 psi 


(38) 


When Poisson’s ratio v is 0.3, the unrelaxed shear modulus G u is 
calculated to be 

E 

- zcrh " 2 - 5 * 106 psl (59) 

The reciprocal of the relaxed shear modulus G r , by formula (15), is 


1 3 (l - 2 v) -1^,-6 --1 

— = — - A — = L = 1.8x10 D psi 


G r E r 




(40) 


Thus, 


G 

G 


u 


r 


4.5 


(41) 
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To find the third constant, use the special form of stress-strain 
relations (20) for uniaxial constant compression, namely: 

o = 2G r (l + v) e + 2Gp(l + v) x^£ (4 2) 

and solve for x a : 

T ° - i [W Vrr - ] (43) 

By taking from the curve of figure 5 several pairs of values of £ and 
£ and by using the values of a and G p just given, an average x a 

value of 30 hours may be obtained. The stress cr that causes failure 
can now be determined as a function of the critical time (or lifetime). 

To find the upper critical stress (<J cr ) u = F^/h, use formula (25) with 

k = 4, the buckling value for a long, simply supported plate. Substitu- 
tion of the numerical values yields 

h 2 « 

(tf r ) = 23.5 *10 6 psi (44) 

cr ft ^ 

The ratio of this critical upper value to the critical value corresponding 
to a lifetime of, say, x is on the basis of the preceding theory the 
same as the ratio of G(x)/g u given by equation (30). Figure 6 shows a 

plot of the critical stress against lifetime for four thickness-to-width 
ratios b/h of 20, 30, 45, and 60. 

These purely theoretical lifetime curves can be compared with ex- 
perimental results presented by Mathauser and Deveikis in reference 3, 
which are also shown in figure 6. Even though the quantitative discrep- 
ancy between the predicted and experimental results is in general fairly 
large, the theory predicts more than merely an order of magnitude and 
exhibits the proper dependence on time and on thickness-to-width ratio. 

Because of the scarcity of information on material properties of 
aluminum alloys subjected to compression at high temperatures and in- 
formation on creep lifetime of plates, this first numerical comparison 
of theoretical and experimental results can be regarded as satisfactory 
and promising. 

Performing the comparison on a wider scale did not appear possible, 
because only for the 2024-T3 aluminum- alloy plates at 450° F were, to 
the authors r knowledge, both material constants and creep lifetimes of 
plates in compression determined experimentally. 
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CONCLUDING REMARKS 

In conclusion, it should be emphasized that the present report con- 
tains a highly simplified and idealized analysis, since the influence of 
metallurgical changes, nonlinearities, time- dependent loads, and so forth, 
were not taken into account. 

However, the fact that creep analysis is bound to be a highly in- 
exact endeavor should be noted. The large amount of scatter to be ex- 
pected may be seen, for example, in reference 13, where it is shown that 
varying the temperature ±10° F every 7 minutes produced creep rates six 
times as large as those produced in a steady temperature test at 1800° F. 
Further, a 15-percent stress variation may change a column lifetime by a 
factor of 10. 

Thus it follows that, from a practical point of view, all that can 
be expected from a creep analysis is merely an order of magnitude. 
Precisely for this reason, any attempt to refine the approach by including 
various effects has to be considered with caution. 


Columbia University, 

New York, N. Y., March 17, 1957 
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Figure 4. - Buckling load against critical time 
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